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ABSTRACT 

We present luminosity and surface-brightness distributions of 40 1 1 1 galaxies with 
A-band photometry from the United Kingdom Infrared Telescope (UKIRT) Infrared 
Deep Sky Survey (UKIDSS) Large Area Survey (LAS), Data Release 3 and optical 
photometry from Data Release 5 of the Sloan Digital Sky Survey (SDSS). Various 
features and limitations of the new UKIDSS data are examined, such as a problem 
affecting Petrosian magnitudes of extended sources. Selection limits in A- and r-band 
magnitude, A-band surface brightness and A-band radius arc included explicitly in the 
1/Vtnax estimate of the space density and luminosity function. The bivariate brightness 
distribution in A-band absolute magnitude and surface brightness is presented and 
found to display a clear luminosity-surface brightness correlation that flattens at high 
luminosity and broadens at low luminosity, consistent with similar analyses at optical 
wavelengths. Best fitting Schechter function parameters for the A-band luminosity 
function are found to be M* - 5 log ft = -23.19 ± 0.04, a = -0.81 ± 0.04 and <f>* = 
(0.0166 ± 0.0008)ft 3 Mpc~ 3 , although the Schechter function provides a poor fit to 
the data at high and low luminosity, while the luminosity density in the K band is 
found to be j — (6.305 ± 0.067) x 10 s L Q ft,Mpc~ 3 . However, we caution that there 
are various known sources of incompleteness and uncertainty in our results. Using 
mass-to-light ratios determined from the optical colours we estimate the stellar mass 
function, finding good agreement with previous results. Possible improvements are 
discussed that could be implemented when extending this analysis to the full LAS. 

Key words: surveys - galaxies: fundamental parameters - galaxies: luminosity func- 
tion, mass function - galaxies: statistics - infrared: galaxies. 



1 INTRODUCTION 

It is possible to learn much about a population's history by 
taking a census of the present-day population. With a large 
and increasing number of deep, large-area surveys taking 
place, a census of the low-redshift galaxy population may 
be undertaken. Deep imaging allows many different prop- 
erties to be studied simultaneously. Much can be learned 
about galaxy formation and evolution by investigating these 
properties: how they correlate with each other, and the sub- 
populations that exist. 

The advantages of working at low redshift (z ~ 0.1) 
are as follows: (1) a more complete sample may be stud- 
ied, including galaxies with low luminosity or low surface 
brightness, (2) the more luminous galaxies may be studied 
in more depth, investigating morphology and structure as 
well as luminosity and colour and (3) the evolution and se- 
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lection effects that plague high-redshift surveys are less of a 
problem. 

The advantages (in principle) of building such a census 
on near-infrared (NIR) observations are well known. First, 
mass-to- light (M/L) ratios in the NIR are largely insensitive 
to galaxy or stellar type, certainly much less than in the op- 
tical (Bell & de Jong 2001). This means that the NIR light 
is a good tracer of the total stellar mass in a galaxy. More- 
over, the range of M/L ratios is much smaller in the NIR, so 
uncertainties in the stellar mass are much smaller. Not only 
does this mean that a survey limited in NIR magnitude will 
be approximately limited in 'apparent stellar mass', but also 
that morphological measurements in the NIR, for example 
the Sersic index and the half-light radius, will reflect the 
distribution of stellar mass within the galaxy, whereas such 
measures in the optical will be significantly biased by the 
presence of young stellar populations. The A-band galaxy 
luminosity function (LF) is, for these reasons, a convenient 
quantity for numerical or semi-analytic models to predict: 
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Table 1. Sample sizes of _ft"-band galaxy luminosity functions. 



Paper 


Number of galaxies in sample 


Mobashcr ct al. (1993) 


181 


Glazcbrook ct al. (1995) 


124 


Gardner et al. (1997) 


510 


Szokoly et al. (1998) 


175 


Loveday (2000) 


345 


Kochanek ct al. (2001) 


3878 


Cole et al. (2001) 


5683 


Huang et al. (2003) 


1056 


Bell et al. (2003) 


6282 


Eke ct al. (2005) 


15 644 


Jones ct al. (2006) 


60 869 


This work 


40111 



see e.g. Croton et al. (2006), Bower et al. (2006), De Lucia 
& Blaizot (2007) and Bertone, De Lucia & Thomas (2007). 

A second advantage is that the K -corrections in the 
K band are also relatively independent of galaxy type 
(Mannucci et al. 2001), leading to smaller uncertainties in 
the absolute magnitudes. 

A third advantage is that dust is much less of a problem 
in the NIR than in the optical. This means that, whereas op- 
tical measurements of galaxy properties are affected by dust 
obscuration, and therefore strongly dependent on the incli- 
nation of the galaxy, producing a smoothing of the galaxy 
LF, this is not such a problem in the NIR (Driver et al. 
2007b; Mailer et al. 2009). 

However, the main disadvantage of the NIR (for ground- 
based telescopes) is the sky brightness, which is around 
13.5magarcsec -2 in K for the data used here (Dye et al. 
2006). 

There have been several studies of the low-redshift K- 
band LF using Two-Micron All-Sky Survey (2MASS) imag- 
ing combined with various redshift surveys (Cole et al. 2001; 
Kochanek et al. 2001; Bell et al. 2003; Eke et al. 2005; Jones 
et al. 2006). Other estimates have been made by Mobasher, 
Sharpies & Ellis (1993), Szokoly et al. (1998) and Loveday 
(2000), using optically selected samples, and by Glazebrook 
et al. (1995), Gardner et al. (1997) and Huang et al. (2003), 
using samples selected in the K band. Table 1 shows the 
sample size of these K-band LF estimates. 

An accurate detemination of the local Tf-band LF also 
has great value as a baseline for comparison with studies 
at higher redshift (e.g. Cirasuolo et al. 2007). The princi- 
pal uncertainty remaining at low redshift is connected with 
the low-luminosity end of the LF. Moreover, there has been 
significant discussion about possible low-surface-brightness 
incompleteness in 2MASS (Andreon 2002), which would af- 
fect the low-luminosity end of the LF. 

Here we present the first statistical study of galaxies 
in the UKIRT Infrared Deep Sky Survey (UKIDSS) Large 
Area Survey, while leaving detailed investigation of surface 
brightness completeness and the very faint-end of the LF to 
future work. 

This paper is organized as follows. Details about the 
sample used are found in Section 2. The method of analysing 
the data is described in Section 3. In Section 4, the bivariate 
brightness distribution (BBD) and LF are presented. The 



stellar mass function (SMF) is estimated in Section 5. There 
is a discussion in Section 6, followed by the conclusions in 
Section 7. 

For ease of comparison with previous results, a flat cos- 
mological model with Q.m = 0.3 and Qa = 0.7 is used, 
with Ho — 100ft km s -1 Mpc -1 . AB magnitudes are used for 
Sloan Digital Sky Survey (SDSS) magnitudes and Vega mag- 
nitudes for A"-band quantities. For reference, AB and Vega 
magnitudes are related in the r band by tab. = rv cg a + 0.146 
and in the K band by Kab = Kv cga + 1.900 (Hewett et al. 
2006). 



2 DATA 

All of the galaxies used in this sample are drawn from the 
main galaxy sample of Data Release 5 (DR5) of the SDSS 
(Adelman-McCarthy et al. 2007), from which the optical 
photometry and spectroscopic redshifts used in the analysis 
below are obtained. The imaging survey covers 8000 square 
degrees, yielding a sample of 783 070 target galaxies. As of 
DR5, the spectroscopic sample had covered 5740 deg 2 of the 
imaging area, with high redshift completeness. 

The UKIRT Infrared Deep Sky Survey (UKIDSS) is de- 
fined by Lawrence et al. (2007). UKIDSS uses the UKIRT 
Wide Field Camera (WFCAM; Casali et al. 2007) and a 
photometric system described by Hewett et al. (2006). The 
pipeline processing and science archive are described by Ir- 
win et al. (in preparation) and Hambly et al. (2008), respec- 
tively. 

UKIDSS consists of five surveys at a variety of depths 
and areas and using various combinations of the ZYJHK 
filters. 7f-band data for this paper are taken from the 
UKIDSS LAS, which is contained within the field of SDSS. 
UKIDSS Data Release 3 (DR3; Warren et al., in prepara- 
tion), which we have used here, was released in 2007 Decem- 
ber, with the LAS containing coverage in YJHK, including 
1189 deg 2 of coverage in K to a 5<r depth of 18.2 mag, in 
directions towards both the north and south Galactic poles 
(NGP and SGP, respectively). 

2.1 Area covered by the UKIDSS-SDSS sample 

The volume, and hence the area, must be well estimated for 
the normalization of the LF and related quantities. Given 
the complex geometry of the overlap region between the two 
surveys, we estimate the area by dividing the number of 
galaxies in the matched sample by the number density of 
sources (per square degree) in the SDSS sample. 

A substantial fraction of the region of overlap between 
the two samples was surveyed in the SDSS Early Data Re- 
lease, so we use the limit of r < 17.6 (de-reddened Petrosian 
magnitude) rather than r < 17.77 used in later versions of 
the SDSS main galaxy selection algorithm. This corresponds 
to 635 320 target galaxies over the whole area. 

The number density is estimated using the total area 
of the SDSS imaging survey (8000 deg 2 ) and the number of 
galaxies targeted for spectroscopy in the SDSS main galaxy 
sample. But, first the size of this sample must be corrected 
for those objects included that are not galaxies. In SDSS 
DR5, the number of target galaxies with good spectroscopic 
redshifts (zConf > 0.8) and with de-reddened Petrosian 
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Figure 1. Sky coverage of the sample, showing the principal 
regions used for the jackknife samples (see Section 3.1), each of 
which is further subdivided into strips in RA (four in the NGP4 
region and five in the others) giving 24 jackknife regions in total, 
each containing approximately the same number of galaxies. In 
the regions with coverage, a darker shade of grey corresponds to 
a higher number density of galaxies. 



Figure 2. A comparison between the catalogue i^-band Pet- 
rosian magnitudes and the magnitudes estimated by interpolat- 
ing between the fixed aperture fluxes, for the 36 660 galaxies in 
the final sample that are not flagged as dcblended. Contours show 
the density of sources on a logarithmic scale, with sources shown 
as points where the density is low. The mean offset is set to be 
zero and the standard deviation is 0.018 mag. 



magnitudes brighter than r = 17.6 is 391 052. Of these, 
384 617 are spectroscopically classified as galaxies. This sug- 
gests that around (391 052-384 617)/391 052 = 1.6 per cent 
of the target galaxies are not in fact galaxies. Taking this 
into account gives a corrected target sample size of 624 865, 
giving a source density of 78.11 galaxies per square degree. 

The two data sets were matched using the WFCAM 
Science Archive (WSA). 1 The initial sample was found by 
selecting all closest matches within 2 arcsec with good, non- 
zero spectroscopic redshifts (SDSS zConf > 0.8), classified 
spectroscopically as galaxies and having no major quality 
control issues flagged in K (UKIDSS kppErrBits < 256), 
yielding a sample of 49 255 galaxies. The sky coverage of 
this matched sample is shown in Fig. 1. 

This number will be affected by failed detections, which 
may introduce some bias into the sample. One type of failed 
detection is when redshifts have not been obtained, due to 
lack of coverage, failed redshifts or lack of available fibres to 
measure the spectra ('fibre collisions'). These are assumed to 
introduce no bias into the sample, although it has been noted 
(Blanton et al. 2003b, 2005) that the SDSS fiber collisions 
lead to slight incompleteness at high-density regions, which 
may introduce a small bias against the type of galaxy found 
in such environments. 

Another kind of failed detection is when there is a 
problem with the UKIDSS catalogue extraction. During the 
course of this work, a bug in the source extraction was dis- 
covered affecting the deblending algorithm. It was found 
that, for Petrosian magnitudes m Y, H and K, the de- 
blender, when invoked, was making the source significantly 
brighter. (The J-band data are micro-stepped, unlike Y, H 
and K, which may explain why this problem is not seen 
for J-band Petrosian magnitudes.) This bug is in the pro- 
cess of being fixed; the reader is referred to the WSA web 



site for updates. 2 However, in the present analysis, for the 
4835 sources flagged as deblended (almost 10 per cent of 
the matched sample), we have estimated the flux within the 
Petrosian aperture by quadratic interpolation between the 
fixed-aperture fluxes, taking advantage of the fact that the 
Petrosian radii, rp, have not been seriously affected by the 
bug (note that 'even substantial errors in rp cause only small 
errors in the Petrosian flux'; Blanton et al. 2001). By com- 
paring the catalogue Jf-band Petrosian magnitudes for non- 
deblended sources with those estimated by this method, we 
find a scatter of 0.018 mag, which suggests that the method 
is sufficiently accurate to provide Petrosian magnitudes for 
those sources affected by deblending. This is shown in Fig. 
2. 

Of the 49 255 sources in the matched sample, 46 sources 
are flagged as having bad pixel(s) in the default aperture. 
These have been removed, leaving 49 209 in the sample. 

Fig. 3 shows the r-band absolute magnitude of galaxies 
in the whole sample, and of sources excluded because of 
problems with the imaging. 

A small number of the remaining galaxies have very 
large magnitude errors, greater or much greater than 
0.15 mag. So, one final cut is to remove sources with large 
uncertainty in magnitude, in order to restrict the systematic 
errors in our results. Of the 49 209 galaxies remaining, those 
with magnitude errors greater than 0.15 mag in r or K are 
removed (69 in K, 811 in r and 2 in both) leaving 48 327 
galaxies. Given the small number affected by this cut, and 
from Fig. 3, it is assumed that any bias induced by this cut 
will be negligible. 

By estimating the area in this way, the assumption is 
that all SDSS target galaxies would be detected in the LAS, 
if that part of the sky has been surveyed. If this is not the 
case, it will have two effects: (1) particular types of galaxies 



1 http://surveys.roe.ac.uk/wsa/ 



2 http : //surveys . roe . ac . uk/wsa/knownlssues . html 
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Figure 3. r-band absolute magnitude of the sources in the entire 
sample ('All', 49 255), the final matched sample ('Good', 48 327), 
those excluded due to poor A"-band imaging, i.e. bad pixels ('kp- 
pErrBits', 46) and those excluded due to uncertainties in the Pet- 
rosian magnitudes ('Magnitude errors', 882). 

will be underrepresented in the sample (those within the 
SDSS completeness limits but outside the limits for the LAS) 
and (2) the overall normalization will be too high, as the area 
and hence the volume probed will be underestimated. 

Fig. 4 shows r-band Petrosian magnitude and the r — K 
Petrosian colour (note that the apertures are not the same 
in r and K so this is not a true colour) for the sources in 
the matched sample. From the figure it can be seen that (1) 
there are likely to be very few sources at all lying within the 
SDSS flux limit but outside the if-band limit, and (2) many 
of these sources are detected anyway, since there are many 
sources detected fainter than the (nominal) A'-band limit. 
This suggests that the effect on the overall normalization 
will be negligible, although the colour-dependent bias will 
be considered later. 

The effective area can now be given as 48 327/78.11 = 
618.72 deg 2 . 

Calculating the area in this way takes into account any 
other random (not bias-inducing) sources of incompleteness 
that have not been considered here explicitly. 

Note that the final sample will be smaller yet, due to 
limits imposed on magnitudes, radius and surface bright- 
ness, but this is related to the redshift limits (which affect 
the volume probed) rather than the area covered. 

2.2 Apparent magnitudes 

Petrosian (1976) apparent magnitudes are used here (for 
comparison with isophotal magnitudes, see Blanton et al. 
2001) and Galactic extinction corrections are used through- 
out, based on the extinction maps of Schlegel, Finkbeiner & 
Davis (1998). 

Fig. 5 shows the number counts as a function of K- 
band Petrosian magnitude for the whole UKIDSS LAS DR3. 
Comparison with previous results shows evidence of incom- 
pleteness in the galaxy number counts fainter than K = 17. 
We choose to limit the sample to K < 16. We use a faint 
magnitude limit of 17.6 in r. Corrections are not applied to 



Figure 4. r-band Petrosian magnitude and r — K Petrosian 
colour for 48 327 galaxies in the matched sample, before imposing 
a cut at K < 16. Objects lying below the solid (red) line have 
K > 16 and will not be included in the final analysis. Contours 
are spaced on a logarithmic scale, with sources shown as points 
where the density is low. 
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Figure 5. i\-band number counts for the full DR3 LAS sample 
(requiring a seamless sample, with no quality control error bits 
flagged), for the final LF sample, and for various samples from 
the literature. The blue crosses show the number counts for the 
whole sample (15 079199 objects), while the red plus-signs are 
for sources classified as galaxies (kClass = 1; 9 962 258 objects). 
Poisson error bars are shown. It may be seen that the counts for 
the LF sample lie always below the counts for the whole UKIDSS 
sample and have a different slope, suggesting that, at all magni- 
tudes, the sample is limited by other factors in addition to the 
A'-band magnitude. 

compensate for the effect of seeing on the Petrosian magni- 
tudes, or for the different fraction of the galaxy's flux recov- 
ered by the Petrosian magnitudes for different galaxy types 
(99 per cent for an exponential profile and 82 per cent, or 
+0.22 mag, for a de Vaucouleurs profile; Blanton et al. 2001). 

Bright limits are applied to the SDSS fibre magnitudes, 
since the spectroscopic sample requires these to be fainter 
than 15 in g and fainter than 14.5 in i. These, however, 
have no effect on our results since the strongest constraint 
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on the low-redshift visibility of each galaxy is generally the 
large radius limit. No bright limit has been applied in the 
if -band; adding a limit as strong as K > 13 has a barely 
noticeable effect. 



2.3 Radii 

Two radii are used below: the Petrosian radius and, for cal- 
culating the surface brightness, the elliptical half-light (semi- 
major axis) radius. 

Given the sky brightness in the if -band and the depth 
of the LAS, the sky subtraction algorithm has to operate 
with a relatively small pixel size of around 24 arcsec. This 
is set in the UKIDSS pipeline as the maximum diameter 
for the (circular) aperture over which the galaxy's flux is 
integrated. A diameter of 24 arcsec corresponds to a Pet- 
rosian radius of 6 arcsec, so this means that any source with 
a true Petrosian radius greater than 6 arcsec will have its 
radius clipped at this value, and the quantity labelled as the 
Petrosian magnitude will in fact be an underestimate of the 
Petrosian flux. 

To avoid complication, a large radius limit has been set, 
requiring the Petrosian radius to be less than 6 arcsec. It is 
worth noting that a significant fraction of the sample (7513 
out of 48 327) have their Petrosian radii clipped at 6 arcsec. 
This limit is taken into account when estimating the contri- 
bution of each galaxy to the space density. However, when 
the large-radius limit is combined with the maximum red- 
shift limit of z = 0.3 (see Section 2.5), it will have the effect 
of excluding all galaxies above a certain physical size. With 
the angular diameter distance, Da{z), given in ft -1 Mpc, the 
Petrosian radius of the largest galaxy that could be visible 
in the sample is 



R = 



10007r(6arcscc)D A (0.3) 
180 x 3600 



= 18.7h -1 kpc. 



(1) 



This will lead to an incompleteness of very large galaxies in 
the final results. 

No small radius limit has been set, as the problem of 
misclassifying galaxies as stars is thought to be negligible in 
SDSS (Blanton et al. 2001), and we do not require indepen- 
dent classification in the K band. 

The scale-size of a galaxy is conventionally measured 
using the radius enclosing half of the total light, known as 
the half-light or effective radius and denoted by r e . In order 
to correct for inclination, the half-light radius we would like 
is the semi-major axis of the elliptical aperture, of the same 
cllipticity and position angle as the galaxy, containing half 
the total flux of the galaxy. This is estimated (Cross et al., in 
preparation) using the Petrosian flux, the 13 circular aper- 
ture fluxes, the ellipticity and the seeing, all made available 
in the WSA (the pipeline does not measure the half-light 
radius). 

2.4 Surface brightness 

The half-light, or effective, surface brightness is estimated 
from the half-light radius by (Blanton et al. 2001) 

2 



fx e — m + 2.5 log 2irr c 



(2) 




band effective surface brightness / (mag arcsec 2 ) 



Figure 6. Number counts as a function of A"-band effective 
surface brightness. The limits in magnitude and radius have been 
imposed on the matched sample, leaving 40 330 sources. Poisson 
errors arc shown. 



al., in preparation) suggest that, for a de Vaucouleurs 
profile, galaxies with surface brightness fainter than 
19.5 mag arcsec -2 are likely to have their fluxes and sizes 
underestimated. For well-defined sample limits, we should 
impose a cut in surface brightness at that value. However, 
when investigating the space density of galaxies with high 
surface brightness, this limit can be safely ignored, since 
for the vast majority of the sample, the faint magnitude 
limit in r provides a stronger constraint on the visibility of 
the galaxy. So in order to include at least some low-surface 
brightness galaxies in the analysis we set a limit in the K- 
band of /i c < 21 mag arcsec -2 . 

The combined limits in if -band Petrosian magnitude 
and Petrosian radius will impose a limit on the effective 
surface brightness for each galaxy. For Sersic (1968) indices n 
between 1 and 4, the Petrosian radius is approximately twice 
the effective radius (Graham et al. 2005), so the faintest 
effective surface brightness will be given by 



He ^ 16 + 2.51og27r(6arcsec/2) 
= 20.38 mag arcsec - 



(3) 
(4) 



where m is the Petrosian magnitude. 

Tests on the UKIDSS source extraction (Cross et 



using the limits in Petrosian magnitude and radius described 
above. Fainter than this, there is a sharp decrease in the 
number counts as a function of surface brightness. We can 
therefore expect to find significant incompleteness at low 
surface brightness. 

Fig. 6 shows the number counts as a function of if -band 
effective surface brightness, after the limits in magnitude 
and radius have been applied. 

The SDSS main galaxy sample has a limit of [i e ^ 
24.5 mag arcsec -2 (Strauss et al. 2002). This limit is taken 
into account here, although it has a negligible effect on 
our results. Very few galaxies in our sample have /i c , r > 
23 mag arcsec -2 , so we assume the limit in SDSS surface 
brightness adds no further incompleteness to the sample, 
once the magnitude limits in r and K and the surface bright- 
ness limit in if have been considered. 
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z log counts 

Figure 7. Redshift and i\-band absolute magnitude distribu- 
tion of the sample (contours and points), with histogram of red- 
shift distribution and (logarithmic) histogram of if -band absolute 
magnitude distribution. For reference, the absolute magnitude as 
a function of redshift corresponding to a source at the if-band 
faint magnitude limit, with typical K- and evolution-corrections 
and neglecting the r-band limit, is shown by the blue dashed 
curve. It can be seen that relatively few galaxies are observed 
near the if-band magnitude limit; this is because of the r-band 
magnitude limit. Redshift distributions for a subdivided sample 
are shown in Fig. 16. 



2.5 Redshifts 

Galactocentric velocity corrections are applied (Loveday 
2000), which typically change the redshifts such that each 
galaxy in this sample is 0.005 ± 0.01 mag fainter, with some 
low-redshift galaxies (z > 0.01) changed by almost 0.1 mag. 

Fig. 7 shows the distribution of the sample in redshift 
and if -band absolute magnitude, excluding those sources 
that lie outside the limits in (apparent) magnitude, radius, 
surface brightness and redshift used in the analysis below. 
The presence of large-scale structure can be seen. 

A high redshift limit of z < 0.3 is imposed to limit the 
effect of K- and evolution-corrections (see below). 

A low redshift limit of z > 0.01 is chosen to limit the 
effect of peculiar velocities, which are not taken into account. 
Note that this limit would need to be relaxed in order to 
sample galaxies with very low luminosity. 

Table 2 shows the various limits on the sample. 



2.6 K- and evolution-corrections 

isT-corrections (to 2 = 0) are estimated using kcorrect 
v4_l_4 (Blanton & Roweis 2007). The five optical Petrosian 
magnitudes (ugriz), before the evolution corrections have 
been applied (see below), are used to fit galaxy templates to 
each galaxy, from which the ./^-correction is derived in r and 
K. WFCAM filter files have been generated from Hewett 
et al. (2006). Due to the inconsistent Petrosian apertures 
between UKIDSS and SDSS bands, we do not have good 
optical-NIR galaxy colours and so we have not been able 
to use the full set of bands (ugrizYJHK) for the template 
fitting. 

When covering a significant range in redshift, it is im- 



Table 2. Limits set on observed quantities, used to define the 
sample and to estimate the contribution of each galaxy to the 
space density. 



Quantity 



Minimum 



Maximum 



K Petrosian magnitude 
r Petrosian magnitude 
g fiber magnitude 
i fiber magnitude 
K Petrosian radius 
Me, if 



15 mag 
14.5 mag 



0.01 



16 mag 
17.6 mag 



6 arcsec 
21 mag arcsec -2 
24.5 mag arcsec -2 
0.3 



portant to include evolution corrections. This is because the 
bright end of the LF will consist mainly of galaxies at high 
redshift, while the faint end will be made up from galaxies at 
low redshift. A failure to include evolution corrections leads 
to a distortion in the shape of the luminosity function. 

Evolution corrections are applied using a simple E(z) = 
Qz, where Q in ugriz is taken from Blanton et al. (2003a) 
and we take Q = 1 in K (consistent with stellar population 
synthesis models, Blanton et al. 2003a, Section 6.2), i.e. Q = 
(4.22,2.04,1.62,1.61,0.76,1.0) for ugrizK. Note that these 
evolution corrections are very simplistic, since different kinds 
of galaxies evolve in different ways, but the correction is 
small: at z = 0.3, the K -band evolution correction is Qz = 
0.3 mag. 

Absolute magnitudes are given by 



M = m- DM{z) - K(z) + E{z) 



(5) 



where DM(z) is the distance modulus and K(z) is the K- 
correction. Absolute surface brightness, in mag arcsec -2 , is 
given by 



Mabs = At - 101og(l +z) - K(z) + E(z) 



(6) 



Note that in the K band, for all galaxies in the sample, the 
absolute surface brightness is not more than 0.5 mag brighter 
than the apparent surface brightness. K- and evolution- 
corrections are not applied when estimating the physical 
radius, which is given in /i -1 kpc as 



R 



10007r6>P A 
180 x 3600 



where 6 is the angular size in arcsec and Da is the angular 
diameter distance in ft -1 Mpc. 

Fig. 8 shows the K- and evolution-corrections in the K 
band and r band, also showing the if-band corrections used 
by Bell et al. (2003). Their ^-corrections are stronger than 
those used here, but their value of Q = 0.8 is weaker than 
our value. These two largely cancel each other out, with our 
K(z) — E(z) being approximately the same as theirs for low 
redshift. 



2.7 Dust 

Driver et al. (2007b) have shown that dust has a consid- 
erable effect on the shape of the optical LF, and caution 
that the effect could still be significant in the NIR. To cor- 
rect for this inclination-dependent dust attenuation would 
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Figure 8. K- (shaded regions) and evolution- (solid lines) cor- 
rections in the r band (upper panel) and in the K band (lower 
panel). Also shown for the K band are the mean i^-correction 
(dashed line, K(z) = — 2.1z) and evolution-correction (dotted 
line, E(z) = 0.8z) from Bell et al. (2003). (These are simple ap- 
proximations, so close agreement is not to be expected.) 

ideally require good bulge-to-disc decomposition, which is 
beyond the scope of this work. We choose here to present 
results that are not corrected for dust, thus representing the 
amount of K-ha,nA light that escapes from the galaxy, rather 
than the amount of A"-band light emitted by the stars in the 
galaxy (some of which will be absorbed by dust). However, in 
estimating the absolute magnitude of a galaxy we have im- 
plicitly assumed isotropic emission of light, which is not the 
case: a disc galaxy will appear fainter when viewed edge-on 
than when viewed face-on, leading to a corresponding under- 
or over-estimate, respectively, of the galaxy's total (attenu- 
ated) luminosity. This will lead to a blurring of the LF for 
disc galaxies, but less so than at optical wavelengths. 



3 METHOD 

In a survey limited by more than one measured quantity, 
e.g. two magnitudes, or magnitude and surface brightness, 
there are two ways to proceed. One is explicitly to take into 
account all of the known selection effects. The other is to set 
strong limits in one or more quantities so that the selection 
effects in the remaining quantities are negligible. Either of 
these ensures that completeness is taken into account and 
that the contribution of each observed galaxy to the total 
space density is estimated accurately (by 1/V max ). 

There are three traditional methods for estimating the 
LF. If the functional form of the LF is known, the best 
method to use is the STY (Sandage, Tammann & Yahil 
1979) maximum likelihood method. But in this work the 
space density is a four-dimensional distribution and the 
functional form is unknown, so the STY method cannot 
be used. The other two methods have no prior assumptions 
about the functional form. These are the stepwise maximum 
likelihood (SWML) method (Efstathiou, Ellis & Peterson 
1988), a version of the STY method where the functional 
form is described by a value in each bin, and the 1/Vmax 
method (Schmidt 1968), which is a simple estimator of the 



value of the space density in each bin. The SWML method 
works by factoring out the normalization at each redshift 
and matching the shape of the LF across different redshifts. 
It is therefore less prone to clustering effects than the 1 /y max 
method, which assumes a homogeneous distribution, but it 
(like the STY method) does not give a value for the overall 
normalization, which must be determined by other means. 
However, the bins in the SWML method are strongly corre- 
lated, so any problem in the SWML estimation is likely to 
affect the entire luminosity function. The 1/Vinax estimator 
is more robust against systemic failure, but requires a large 
survey volume in order to overcome the effects of large-scale 
structure. 

In this work, due to the complexity of the four- 
dimensional space, we have found the SWML method to 
be prone to fail, whereas the 1/V max method has proved to 
be more stable. Large-scale structure has not been taken 
into account, but a large volume has been probed at bright 
luminosities and faint luminosities are already affected by 
incompleteness at low surface brightness. 

The volume within which each galaxy would be visible, 
Vmax, is estimated for each galaxy by considering the tightest 
constraints on maximum and minimum redshift provided 
by the limits in z, K, r, 7f-band fi B and K-band Petrosian 
radius. Each galaxy is then given a weight u>i = l/Knax.i. 
The space density <f) of galaxies of a certain type (e.g., binned 
in absolute magnitude and/or surface brightness) is given by 
<f> — Wi for galaxies i of that type. 

The effect of the 1/V max factors on the features of the 
AT-band LF may be appreciated by comparing the shape of 
the absolute magnitude histogram in Fig. 7 with the shape of 
the LF in Fig. 14. For galaxies with higher luminosity, V max 
is considerably greater, so the effect of the weighting on the 
shape of the curve is to raise the faint end relative to the 
bright end. It may also be seen that the 1/V max weighting 
introduces no obviously spurious features into the final LF, 
despite the complexity of the limits on the sample. 

Some regions of the four-dimensional parameter space 
have not been sampled, for example, low-surface brightness 
galaxies, so care needs to be taken when interpreting the 
binned space density. 

3.1 Errors 

Statistical errors for all quantities (except where indicated) 
are estimated through the jackknife method (Lupton 1993). 
The sample area is divided into 24 regions, each containing 
approximately the same number of galaxies; the principal 
subdivisions are shown in Fig. 1. The space density is calcu- 
lated once for the whole sample and then a further 24 times, 
each time omitting one of the 24 regions. The variance of 
each value of (j> is then calculated using: 

Var(0) = — Y^iti-W? (8) 

i 

where (f>i is the ith jackknife resampling of the data and n 
is the number of jackknife resamplings (24 in this case). A 
bias correction (Lupton 1993) is applied to the space density, 
giving a new value of 

f = 0+(n-l)(0-£?) (9) 
where <j) is the original estimate for the whole sample. 
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Uncertainties and bias corrections for the functional fits 



and integrated quantities (for example, luminosity density) 
are estimated using the same method. 

It would be preferable, but more difficult, to divide the 
sample into regions of equal area rather than equal num- 
bers of galaxies. Given the way in which the effective area 
is calculated, these subdivisions cannot say anything about 
the overall normalization of the space density; for example, 
whether the south Galactic pole has a lower density of galax- 
ies than the north Galactic pole. 

Uncertainties in the magnitudes and other galaxy prop- 
erties are not included in the analysis. This is likely to have 
only a small effect on most of our results, but the magnitude 
errors will cause us slightly to overestimate the space den- 




sity of galaxies at the bright end of the luminosity function, 0.0195 0.0200 0.0205 

where the LF is very steep (Jones et al. 2006). 



3.2 Functional fits 

Functional fits for the BBD, LFs and SMF are found using 
the mpfit routines written by Craig Markwardt. 3 

For the LFs and the SMF, the Schechter function is 
used. Expressed in magnitudes, this takes the form 

</>(M) = 0.41nlO<A0 ' 4(A/ *- M)(Q+1) e- lo0 - 4<M *~ M) (10) 

where M* is the characteristic absolute magnitude, a is the 
faint end slope and <f>* is the normalization. 

For the BBD, the Choloniewski (1985) function is used 
(Driver et al. 2005): 



0(M lMe ) = ^!Vi0°' 4(M *- M)(Q+1) < 



n 0.4(M* —Ml 



x exp 



ti e -{J.* e -p(M-M*) 



(11) 



where M* , a and cf>* are the usual Schechter function pa- 
rameters, fil is the mean surface brightness at M* , <r Me is the 
standard deviation in surface brightness and j3 is the slope 
of the relationship between absolute magnitude and mean 
surface brightness. Note that integrating the Choloniewski 
function over surface brightness gives a Schechter function. 

The range of points over which the functions are fit is 
restricted in order to avoid known regions of incompleteness: 
only data brighter than —20 in Mk — 51og/i, —18 in M r — 
5 log ft and 19 in surface brightness, or more massive than 
9.5 in log(/i -2 M0), are used to generate the functional fits. 



3.3 Testing the method with simulated data 

We test the V ma x estimator and fitting routines using four 
simulated samples drawn from a known Choloniewski func- 
tion, with Gaussian r — K colours (mean and standard devi- 
ation derived from the observed sample, but with no depen- 
dence on luminosity), and subject to the same observational 
limits of the data sample, except with r < 17.77 rather than 
r < 17.6. Each sample contains between 42 632 and 42 885 
galaxies. 



http: //cow. physics .wise . edu/~craigm/idl/f itting.html 



Figure 9. Schechter function parameters for the simulated sam- 
ples. Plus symbols show the input parameters, M* — 51og/i = 
—23, a = —1 and <f>* = 0.02h 3 Mpc -3 , while contours show the la 
and 2<t error contours on the recovered parameters for the four dif- 
ferent simulated samples. These uncertainties are estimated from 
24 jackknife estimations of the Schechter function. 



Fig. 9 shows the recovered Schechter function param- 
eters from the simulated samples. There is no obvious sys- 
tematic bias. The recovered uncertainties, shown by the con- 
tours, give an approximate measure of the closeness of the 
recovered parameters to the input Schechter function, al- 
though they tend to underestimate the true error. 

The luminosity density, in solar units, may be calculated 
from the Schechter function as 



J = 



L(f>(L) dL = <^*10 - 4(Me - M * ) r(a + 2) 



(12) 



where Mq is the absolute magnitude of the sun (taken to be 
3.32 in the Tf-band, Bell et al. 2003), or from the weights of 
each galaxy by 



4(M -M;) 



(13) 



For the simulated samples the luminosity density was esti- 
mated in these two ways and compared to the luminosity 
density from the input Schechter function. By summing the 
galaxy weights, the recovered luminosity density typically 
underestimated the input luminosity density by around 1 
per cent, whereas there was no obvious bias from integra- 
tion of the recovered Schechter function. 

Fig. 10 shows the input Choloniewski function and 
the recovered BBD for one of the simulated samples. The 
Choloniewski function gives a good fit to the recovered BBD, 
but the recovered BBD is itself biased with respect to the 
input Choloniewski function, for all four simulated samples, 
most noticeably towards higher surface brightness, by ap- 
proximately 0.15 mag arcsec~ 2 . The BBD presented below 
should therefore be considered only approximately correct. 
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Figure 10. Input Choloniewski function for the sim- 
ulations (red dotted contours) and the recovered BBD 
from one of the simulated samples (shaded regions and 
black solid contours). The space density, 0, is in units 
of h s Mpc -3 mag - 1 (mag arcsec -2 ) -1 . A Choloniewski func- 
tion was fit to the recovered BBD (not shown). The in- 
put (recovered) Choloniewski function parameters are M* — 
5 log h = -23(-23.03 ± 0.01) mag, a = -l(-0.96 ± 0.01), 
<j>* = 0.02(0.0224 ± 0.0003)/i 3 Mpc" 3 , /j,* = 17.5(17.363 ± 
0.005) mag arcsec -2 , <r Mo = 0.6(0.580 ± 0.003) mag arcsec -2 and 
(3 = 0.3(0.287 ± 0.003). Very similar results are obtained for the 
other three simulated samples. 



4 BIVARIATE BRIGHTNESS DISTRIBUTION 
AND LUMINOSITY FUNCTIONS 

In this section the bivariate brightness distribution and lu- 
minosity functions are calculated for the whole sample and 
for subdivisions of the data. 



4.1 Bivariate brightness distribution 

Fig. 11 shows the BBD in Tf-band absolute magnitude and 
absolute effective surface brightness. The value of the space 
density at any point on the BBD is estimated assuming the 
visibility of the full range of galaxy types that exist with that 
absolute magnitude and surface brightness. This may in fact 
not be the case, given the additional limits in faint r-band 
apparent magnitude and large T^-band Petrosian radius. 

Figs. 12 and 13 show the effect of the combined r- and 
_K"-band flux limits on the completeness, as a function of K- 
and r-band absolute magnitude respectively. 

From Fig. 4 it can be seen that near the r-band flux limit 
blue galaxies cannot be seen, and that near the K-hanA flux 
limit red galaxies cannot be seen. For the red galaxies that 
would not be seen at near the if-band flux limit, we consider 
galaxies at redshift z with r = 17.6 and various r—K colours, 
neglecting K- and evolution-corrections. These sources will 
have K-h&nA absolute magnitude 

M K =K- DM(z) =r-(r-K)- DM(z) . (14) 

Fig. 12 shows this relation between Mk and r — K (with 
r = 17.6) for z — 0.01 and 0.02, corresponding respectively 
to the minimum redshift considered here and the redshift at 
which the survey volume is just over 10 4 /i- 3 Mpc 3 . It can 
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Figure 11. BBD for the full sample in i^-band absolute 
magnitude and absolute effective surface brightness. Shaded re- 
gions and solid black contours show the space density, (f>, as in 
Fig. 10. The best-fitting Choloniewski function, estimated using 
Mk — 51og/i < —20 and /i c ,a.bs < 19, is shown by the red dotted 
contours. Parameters of the fit are M* — 51og/i = —22.96 mag, 
a = -0.38, <j>* = 0.0201/i 3 Mpc" 3 , fi* abs = 17.36 mag arcsec -2 , 
<r Mc aba = 0.672 mag arcsec -2 and j3 = 0.188. 



og 




-18 -20 -22 -24 -26 

M K -5log h 



Figure 12. Space density of galaxies, in units of 
h 3 Mpc -3 mag -2 , as a function of A"-band absolute magnitude 
and rest-frame r — K Petrosian colour. Note that this is not a 
true colour since the apertures differ between the bands. The 
straight lines show the position on the plot of hypothetical galax- 
ies at the r-band flux limit, with various r — K colours, situated 
at z = 0.01 (upper-left line, solid, red) or at z = 0.02 (lower-right 
line, dashed, green), z = 0.02 corresponds to a survey volume of 
1.17 X 10 4 h -3 Mpc 3 . 

be seen that red galaxies will not be seen at faint A'-band 
luminosity, and we therefore anticipate some incompleteness 
at M K - 51og/i > -19. 

A similar relation may be derived for the blue galaxies 
that would not be seen near the isT-band flux limit. Consid- 
ering a galaxy at redshift z with K = 16 (the K-band flux 
limit), the r-band absolute magnitude will be 

M r = r- DM(z) =K+(r-K)- DM(z) . (15) 
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Figure 13. As Fig. 12 but for r-band absolute magnitude. The 
straight lines correspond to hypothetical galaxies at the A'-band 
flux limit, with various r — K colours, situated at z = 0.01 (lower- 
left line, solid, red) or at z = 0.02 (upper-right line, dashed, 
green). 

Fig. 13 shows this relation between M r and r — K (with 
A — 16), from which it may be seen that blue galaxies will 
not be seen at faint r-band luminosity. Significant incom- 
pleteness is to be expected for M r — 5 log h > —17.5. 

For the large radius limit, considering a source with a 
certain if -band magnitude and surface brightness, and as- 
suming the Petrosian radius rp is twice the effective radius, 
Equation (2) can be written as 

Mc ~A + 2.51og27r(^) 2 . (16) 

For A=16 and fi e = 19.5 mag arcsec -2 , this corresponds to 
a limit in Petrosian radius of 4.0 arcsec, considerably less 
than the 6 arcsec limit intrinsic to the data; only for surface 
brightness fainter than 20.38 mag arcsec" 2 does the 6 arcsec 
large radius limit dominate. 

Fig. 11 also shows the best-fitting Choloniewski fit to 
the BBD. It can be seen that the function provides a poor fit 
to the data, being unable to fit simultaneously the decline 
at high surface brightness and high luminosity, the broaden- 
ing of the surface brightness distribution at faint luminosity, 
or the slope of the luminosity-surface brightness relation, 
which flattens at high luminosity. These features have also 
been seen in optical determinations of the BBD (e.g., Driver 
1999; Blanton et al. 2001; Cross & Driver 2002; Driver et al. 
2005) 

4.2 A'-band LF and luminosity density 

Fig. 14 shows the A'-band LF for the whole sample. The 
parameters of the best-fitting Schechter function correlate 
strongly, with corr(M*,a) = 0.92, corr(M*,<^)*) = 0.97 and 
corr(a,0*) = 0.92. 

The most significant deviation from the published A'- 
band luminosity functions is at the bright end, where our 
LF has a very steep decline compared with those of Bell 
et al. (2003), Eke et al. (2005) and Jones et al. (2006). There 
are several possible explanations for this. For example, the 
following explanations are given. 



(i) Differences in the evolution corrections used (0.3 mag 
at z = 0.3 for our value of Q = 1), affecting the high- 
luminosity galaxies (cf. Fig. 7). Jones et al. (2006) applied 
no evolution corrections, although they have a redshift limit 
of 2 < 0.2. Eke et al. (2005), for z < 0.12, and Bell et al. 
(2003) applied combined A- and evolution-corrections simi- 
lar to our own. There is better agreement between our results 
and those of Kochanek et al. (2001), although they applied 
no evolution corrections. 

(ii) We have used Petrosian magnitudes, which are signif- 
icantly fainter than total magnitudes for galaxies with a high 
Sersic index (Graham et al. 2005), for example, by 0.22 mag 
for a de Vaucouleurs profile (Blanton et al. 2001). Kochanek 
et al. (2001) used isophotal magnitudes with a correction of 
0.20 ± 0.04 mag to estimate total magnitudes, but a larger 
correction may be needed for the most luminous galaxies 
(Blanton et al. 2001), which tend to have a higher Sersic 
index. Bell et al. (2003) used the 2MASS Kron magnitudes 
with a correction of 0.1 mag, Eke et al. (2005) applied a simi- 
lar correction based on the J — As colour, while Jones et al. 
(2006) used the 2MASS total (extrapolated) magnitudes. 
However, this effect is countered to an extent by the effects 
of seeing, which, for a de Vaucouleurs profile, increases the 
fraction of the galaxy's flux recovered by the Petrosian mag- 
nitude when the angular size of the galaxy is small (Blanton 
et al. 2001). This is the case for most of the luminous galax- 
ies in the sample, which are generally observed at higher 
redshift. 

(iii) Unidentified sources of incompleteness or a poorly 
understood selection function, given the non-trivial limits on 
our sample and the dependence on completeness in SDSS. 

(iv) Improved photometry: significant magnitude errors 
at the bright end will lead to an overestimate of the space 
density at high luminosities (Jones et al. 2006). The results 
of Kochanek et al. (2001), Bell et al. (2003), Eke et al. (2005) 
and Jones et al. (2006) are all based on 2MASS magnitudes, 
with much shallower imaging than the UKIDSS LAS. More- 
over, the most luminous galaxies in our sample are all ob- 
served at magnitudes much brighter than the A-band flux 
limit, as can be seen from Fig. 7, so the photometric errors 
are very small. 

Fig. 15 shows a comparison between A-band 2MASS 
Kron magnitudes and UKIDSS Petrosian magnitudes for 
sources in this sample with counterparts in the 2MASS ex- 
tended source catalogue (XSC) (Jarrett et al. 2000). It can 
be seen that a significant number of faint 2MASS sources 
have 2MASS Kron magnitudes that are much brighter than 
UKIDSS Petrosian magnitudes, sometimes by over 0.5 mag. 
If the UKIDSS magnitudes are more accurate, then this sug- 
gests that 2MASS luminosity functions will over-estimate 
the bright end, since high luminosity galaxies are most often 
found near the faint magnitude limit. For a more in-depth 
comparison between UKIDSS and 2MASS photometry, see 
Cross et al. (in preparation). 

(v) Missing sources in high-density environments. As dis- 
cussed in Section 2.1, the positioning of the SDSS spectro- 
scopic fibres leads to a slight incompleteness at high-density 
regions. These regions are rich in high-luminosity galaxies, 
so this may affect the high-luminosity end of the LF. 

(vi) Incompleteness at very large physical radius. As 
shown in Section 2.3, the limits in A-band Petrosian radius 
and redshift combine to set an upper limit to the physical 
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Figure 14. i\-band LF for the whole sample, with a compendium of published results from observations or semi-analytic models. 
Only the filled points are used in the Schechtcr function fit, i.e., Mk — 51og/i < —20; the unfilled points are likely to suffer from some 
incompleteness of low-surface brightness or red, low-luminosity galaxies. Schechter function parameters are M* — 5 log ft = —23. 19 ±0.04, 
a = —0.81 ± 0.04 and cf>* = (0.0166 ± 0.0008)ft 3 Mpc -3 , although, as noted in the main text, the Schechter function provides a poor fit 
to the LF at both high and low luminosity. 



radius of galaxies in the sample. It appears likely that this 
will lead to some incompleteness at high luminosities. 

There is general agreement at the bright end of the (Pet- 
rosian) r-band LF (see below), which is consistent with the 
possible explanations listed above, although it does suggest 
that the effect of missing galaxies is not the most significant 
factor, whether these are large galaxies or galaxies found in 
high-density environments. Further investigation is required 
to resolve this. 

The steeper bright end leads to a fainter value of M* 
compared with previous results, and consequently (given the 
strong correlation between M* and a) a shallower value for 
the faint-end slope, a. It should be clear from Fig. 14 that 
the faint end of the Schechter function does not coincide with 
the faint end of the LF. It should also be noted that the 
1/V m ax method is sensitive to large-scale structure, which 
will be particularly apparent at the faint end of the LF. 

At intermediate magnitudes, Mk — 5 log ft ~ —22, our 
LF is noticeably higher than those from the literature. This 
is likely to be due to an overdensity of galaxies at z < 0.1, as 
seen in Fig. 7. The effect of large-scale structure is discussed 
in Section 4.3. 

The 7f-band luminosity density, with jackknife errors, 
is found to be j = (6. 185 ±0.072) x 10 s L ft Mpc -3 from the 
Schechter function, or j = (6.305 ± 0.067) x 10 8 L Q ftMpc -3 
from the galaxy weights. The true luminosity density is likely 
to be higher than these values since the Schechter function 
provides a poor fit at both the faint end and the bright 



end, and since there are known sources of incompleteness 
in absolute magnitude, colour and surface brightness at the 
faint end. For comparison, measurements of the luminosity 
density from the published results shown in Fig. 14 lie within 
approximately (5.8-7.6) x 10 8 Lq ft Mpc -3 (see Jones et al. 
2006, fig. 15)! 

4.3 Large-scale structure 

In order to show the effect of large-scale structure on our 
results, we have estimated the LF for the five principal sub- 
divisions of the sample (see Fig. 1). 

Fig. 16 shows the K-band LF for the SGP region and 
the four NGP regions, together with the LF for the whole 
sample and the redshift distribution for each sub-sample. 
We note the following. 

(i) There is a large scatter at low luminosities, illustrating 
the limitations of the 1/Vinax method. 

(ii) Each sub-sample is normalized according to the num- 
ber of galaxies rather than the area covered, so a genuine 
over-density at a certain redshift will be artificially compen- 
sated for by an apparent under-density at other redshifts, 
and vice versa. Given the correlation between redshift and 
absolute magnitude (Fig. 7), this means that large-scale 
structure at intermediate redshifts will distort the LF at 
both ends. This may be seen for the SGP region, where 
an under-density at z ~ 0.1, seen at Mk — 5 log ft between 
—22.5 and —23, may have led to an over-estimate of the LF 



12 A. J. Smith, J. Loveday & N. J. G. Cross 




11.0 11.5 12.0 12.5 13.0 13.5 

2MASS K Kron 



Figure 15. Residual if -band magnitudes between 2MASS Kron 
and UKIDSS Petrosian magnitudes for 6898 galaxies matched be- 
tween our sample and the 2MASS XSC, with 2MASS A K roii < 
13.5 mag. Dust reddening corrections have not been applied. Con- 
tours show the density of sources on a logarithmic scale, with 
sources shown as points where the density is low. The mean off- 
set is 0.19 mag, such that 2MASS magnitudes are brighter on 
average, with a standard deviation of 0.29 mag. 
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Figure 16. K-band LF for the NGP and SGP regions, also show- 
ing the LF for the whole sample. Jackknife errors are shown for 
the whole sample and Poisson errors for the sub-samples. The in- 
set shows the redshift distribution of galaxies in each sub-sample, 
with the line colours corresponding to the colours of the symbols 
in the main figure. 



at both the faint end and the bright end, while for the NGP2 
region, an over-density at z < 0.1 {Mk — 5 log h ~ —22.5), 
the Sloan Great Wall (Gott et al. 2005), may have led to an 
under-estimate of the LF at both ends. 

(iii) An over-density affecting the whole sample at a cer- 
tain redshift will be seen at the absolute magnitude cor- 
responding to that redshift (see Fig. 7). There is a large 
over-density at z < 0.1, which may explain why the LF in 
Fig. 14 shows a higher density at Mk — 51og/i between —22 
and —23 compared with the findings of other authors. 
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Figure 17. A'-band LF for galaxies in three different redshift 
ranges, also showing the LF for the whole sample. 



4.4 Variation with redshift 

In order to identify further sources of bias we split the sam- 
ple into three bins in redshift, containing approximately 
equal numbers of galaxies. 

Fig. 17 shows the if -band LF estimated for the low-, 
mid- and high-redshift sub-samples. In the mid- and high- 
redshift sub-samples, there is a gradual downturn of the LF 
towards faint absolute magnitudes. This is to be expected, 
since the faint A"-band absolute magnitude limit at any given 
redshift is strongly dependent on the r — K colour (see Fig. 
12). There is disagreement at the bright end between the 
different slices, with (j> getting progressively higher at higher 
redshift. This trend could be a result of the evolution cor- 
rections being too small, but the trend is still present when 
a correction as strong as Q — 2 is applied. 

Contributing factors could be: (1) an over-simplistic 
form for the evolution corrections, E(z) = Qz, independent 
of galaxy type, (2) a decreasing apparent Sersic index with 
decreasing angular size, as a result of convolution with the 
PSF, causing a greater fraction of a galaxy's flux to be re- 
covered by the Petrosian magnitude when the galaxy is ob- 
served at higher redshift, (3) poorly understood limits to the 
sample, or (4) large-scale structure affecting both the shape 
and normalization of the LF for each redshift slice. 

However, the most significant factor is likely to be the 
limit on the physical size of galaxies included in the sample, 
which arises through a combination of the large angular size 
limit and the high redshift limit (see Section 2.1). In the 
medium- and low-redshift bins, the upper limit to the red- 
shift is lowered, and hence the upper limit to the physical 
Petrosian radius is also lowered. From Equation 1, the phys- 
ical Petrosian radius can be no larger than 9.73ft -1 kpc in 
the medium-redshift bin, and no larger than 6.70ft _1 kpc in 
the low-redshift bin. This will cause significant incomplete- 
ness at the bright end of the LF, since higher-luminosity 
galaxies tend to be larger. 

This test does make it clear that our results are to an 
extent dependent on the redshift limits chosen: choosing a 
lower value than 0.3 for the maximum redshift would have 
given an even steeper bright-end slope for the LF. 



K-band galaxies from UKIDSS LAS 13 



log 



o 3 




• This work 1 /V mQ1 

X Blanton et ol. (2005) SWML, incompleteness corrected 
Blanton et ol. (2003) 



M r -5log h 




Figure 18. r-band LF. Only the filled points are used in the 
Schechter function fit, i.e., M r — 51ogh brighter than —18; the 
unfilled points are likely to suffer from some incompleteness of 
low-surface brightness galaxies or low-luminosity blue galaxies. 
The LF of Blanton et al. (2003a) has been adjusted from the 01 r- 
band at z = 0.1 to the r-band at z = by (1) shifting Qzmag 
fainter in magnitude, with Q = 1.62 and z = 0.1, (2) reducing 
the number density by OAPz dex, with P = 0.18, and (3) making 
the LF 0.22 mag brighter in magnitude, to convert from r to r. 
Schechter function parameters are M* — 5 log h = —20.40 ± 0.04, 
a = -0.91 ± 0.05 and </>* = (0.0195 ± 0.0010)?i 3 Mpc -3 . 



4.5 r-band luminosity function and luminosity 
density 

Fig. 18 shows the r-band LF. While there is excellent agree- 
ment at the bright end with the LF of Blanton et al. (2003a), 
our LF is over-dense at intermediate redshifts, probably due 
to large-scale structure. The deficit of blue low-luminosity 
galaxies, identified in Fig. 13, is clearly evident. 

The r-band luminosity density is found to be j — 
(1.930 ± 0.027) x lO 8 L /iMpc" 3 by extrapolating the 
Schechter function, or j = (1.931 ± 0.019) x 10 s L Q /iMpc" 3 
from the galaxy weights, assuming a solar absolute mag- 
nitude of 4.64 (Blanton & Roweis 2007). Again, the true 
luminosity density is likely to be higher than these values 
given the incompleteness at the faint end. This is some- 
what higher than the r-band z — luminosity density of 
Blanton et al. (2003a), -15.90 + 2.5 log h mag in a Mpc 3 , or 
1.64 x 10 s L Q h Mpc -3 . 

Given that we use the same source of data as Blanton 
et al. (2003a), but over a smaller area, with a more complex 
selection function and with an inferior LF estimator, the r- 
band results we find should not be interpreted as being more 
than a consistency check on our analysis. 



4.6 Subdividing by colour 

The bimodality of the galaxy population has been recog- 
nized by many authors (see Driver et al. 2006; Ball et al. 
2006, and references therein). This may be visualized by 
subdividing the LF or BBD in various ways, for example, 
according to colour, concentration or spectral class. Of these 
properties, we find, following Driver et al. (2006), that the 
u — r core (PSF) colour gives a particularly sharp dichotomy. 



Figure 19. BBD for red galaxies, with (u — r)psF > 2.35. 
The best-fitting Choloniewski function, estimated using Mjf — 
51og/i < —20 and /i e . a b s < 19, is shown by the red dotted 
contours. Parameters of the fit are M* — 51og/i = —22.91 mag, 
a = 0.14, <f>* = 0.0119/1 3 Mpc" 3 , abs = 17.11 mag arcsec -2 , 
<Tu . = 0.584 mag arcsec -2 and (3 = -0.0006. 
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Figure 20. BBD for blue galaxies, with (u — r)psp < 2.35. 
The best-fitting Choloniewski function, estimated using Mk — 
5 log h < —20 and /i c . a b s < 19, is shown by the red dotted 
contours. Parameters of the fit are M* — 51og/i = —22.63 mag, 
a = -0.94, </>* = 0.0115/1 3 Mpc -3 , fi* abs = 17.91 mag arcsec -2 , 
<7u , = 0.851 mag arcsec -2 and f) = 0.419. 



Figs. 19-21 show the K-band BBD and K- and r-band 
LFs split by the SDSS rest-frame u—r PSF colour. The space 
density is estimated by summing the weights of galaxies with 
u — r > 2.35 or u — r < 2.35 for red and blue galaxies 
respectively, with jackknife errors estimated subsequently. 

The BBD for red galaxies, excluding outliers, shows no 
evidence of a correlation between luminosity and surface 
brightness, while the BBD for blue galaxies shows no flatten- 
ing off of the luminosity-surface brightness relation at high 
luminosities, suggesting that this division reflects a property 
of the underlying population. Moreover, the Choloniewski 
function appears to fit the blue BBD much better than it fits 
the BBD for the whole sample. However, we caution that the 
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Figure 21. K- (left-hand panel) and r-band LFs for red and 
blue galaxies, showing the total LF as well. 

lack of red-core galaxies with Mk — 51og/i > —20 could be 
a symptom of the incompleteness identified in Fig. 12, while 
the lack of such galaxies with p e ,abs > 19.5 mag arcsec -2 
could be due to the low-surface brightness limit for de Vau- 
couleurs profile galaxies. 

The LFs show a sharp division, with red-core galax- 
ies more abundant than blue-core galaxies by an order of 
magnitude at high luminosity (and vice versa at low lumi- 
nosity), and with the bright end of the LF around 0.8 mag 
more luminous in the if -band for red-core galaxies. 



Figure 22. Stellar mass function. Only the filled points are 
used in the Schechter function fit, i.e., stellar mass greater than 
lO 9 - 5 h -2 M ; the unfilled points are likely to suffer from some 
incompleteness of low-surface brightness galaxies. Masses calcu- 
lated from the K-band KCORRECT M/L ratios have been in- 
creased by 0.1 dex. Schechter function parameters are found to 
be log{M*h 2 /M e ) = 10.45 ± 0.02, a = -1.03 ± 0.04 and 
0* = (0.0105 ± 0.0006)h 3 Mpc~ 3 . Stellar masses based on other 
IMFs have been reduced for comparison with our results, based 
on the Chabrier IMF: Salpeter IMF (Panter, Heavens & Jimenez 

2004) by 0.3 dex, 'diet' Salpeter (Bell et al. 2003) by 0.15 dex, 
and no conversion has been applied for the Kennicutt (Eke et al. 

2005) or Kroupa (Baldry, Glazebrook & Driver 2008) IMFs. 



5 STELLAR MASS FUNCTION 

Stellar masses are derived from the KCORRECT template fits, 
which are based on the ugriz photometry. Described in more 
detail by Blanton & Roweis (2007), these templates are gen- 
erated from Bruzual & Chariot (2003) stellar population 
synthesis models with a Chabrier (2003) initial mass func- 
tion (IMF). 

The M/L ratio varies less at NIR than at optical wave- 
lengths (Bell & de Jong 2001) so, assuming the uncertainty 
in stellar mass is dominated by the uncertainty in the M/L 
ratio, it makes sense to estimate the stellar mass from the 
NIR absolute magnitude and the M/L ratio at that wave- 
length. Ideally the M/L ratio could be found by fitting a 
template to all available photometry, i.e., ugrizY J H K for 
SDSS and UKIDSS. However, for this analysis, where we 
have good optical colours (with consistent apertures) but 
poor NIR and optical-NIR colours, we find it is best to esti- 
mate the if-band M/L ratio using the optical colours only. 

Fig. 22 shows the SMF, with stellar masses estimated 
from the if -band absolute magnitudes and the A"-band M/L 
ratios from KCORRECT. At the high-mass end our results 
agree well with previously-published SMFs, while at the low- 
mass end the discrepancy could be a result of incompleteness 
or large-scale structure. The underdensity at intermediate 
masses (M — W 9 Ii~ 2 Mq) could be due to inappropriate 
M/L ratios or large-scale structure. 

Fig. 23 shows the SMF calculated using the default mass 
and the various M/L ratios, all given by KCORRECT. We 



found there to be an offset in the masses derived from the K- 
band M/L ratios compared with the masses derived from the 
optical bands. The precise cause of this offset is not known, 
but it has been compensated for by increasing the if -band- 
derived masses by 0.1 dex. The greater uncertainty in blue 
M/L ratios would be expected to stretch the high-mass end 
of the SMF, as seen in the u- and p-bands SMFs. Smaller 
uncertainty in the red or near-infrared M/L ratios may be 
responsible for the disagreement at low masses. However, 
there is some uncertainty in the if-band M/L ratios caused 
by the emission from thermally-pulsating asymptotic giant 
branch (TP-AGB) stars, which are very difficult to model 
(Maraston 1998; Bruzual 2007). 

The stellar mass density is found to be (3.02 ± 0.05) x 
10 s /iMq Mpc -3 by extrapolating the Schechter function or 
(2.99 ± 0.04) x 10 s ft M Q Mpc" 3 from the galaxy weights. 
Due to incompleteness this is likely to be an underestimate, 
and a different IMF could increase this substantially, for 
example, by a factor of 2 for a Salpeter IMF. Baldry et al. 
(2008) find a general consensus that 'the cosmic stellar mass 
density is in the range 4-8 per cent of the baryon density 
(assuming Qb = 0.045)'. In their units, taking h — 0.7, we 
find a stellar mass density of 4.9 ±0.08 or 4.9 ±0.06 per cent 
of the baryon density, respectively, in good agreement with 
previous findings. 
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Figure 23. SMF, using various M/L ratios and the default 
KCORRECT mass, which is derived from the template fit to the 
input (ugriz) absolute magnitudes. Masses calculated from the 
A'-band KCORRECT M/L ratios have been increased by 0.1 dex. 

6 DISCUSSION 

Having already compared our results with those of other 
authors, we now consider some implications of our findings 
along with possible improvements and extensions. 

6.1 Functional fits 

The properties and evolution of the galaxy population 
are conventionally quantified using simple functional (e.g. 
Schechter function) fits to the data. With a large sample 
of low-redshift galaxies, it is possible to test for expected 
bias when such a simple form is assumed for the LF at high 
redshift. 

Fig. 24 shows the Schechter function fits to our _R"-band 
LF, restricting the fit to various portions of the LF. This 
is intended to mimic surveys at higher redshift, where only 
the bright end of the LF is visible. These results suggest 
that, even for a non-evolving LF, with increasing redshift one 
would observe (1) a steeper faint-end slope, (2) a brighter 
characteristic magnitude and (3) a decreasing number den- 
sity at the characteristic magnitude. This arises because of 
systematic deviations from the Schechter function form and 
highlights the danger of relying too strongly on the Schechter 
function fits. 

6.2 Possible improvements and extensions 

The results presented here use just over an eighth of the 
projected final overlap area between the UKIDSS LAS and 
SDSS (4000 deg 2 ). Already the number of galaxies is com- 
parable to the largest A'-band surveys of the low-redshift 
universe, so we now consider how this work could be refined 
for a larger sample. 

The 1/Knax estimator has been shown to be sensitive 
to large-scale structure at low redshift. With a larger sam- 
ple this will be less apparent, but a better estimator of the 
space density is required, e.g., a multivariate version of the 
SWML method or a correction for variation with redshift 
(e.g., Cross et al. 2001). We have found the SWML method 
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Figure 24. A'-band LF, showing Schechter function fits restrict- 
ing the fitting range in absolute magnitudes, intended to mimic 
surveys at higher redshifts. 



difficult to implement over four dimensions, due to the com- 
plexity of the parameter space, but these problems may not 
be insurmountable. 

In this work the uncertainties in the measured quanti- 
ties have been largely ignored and assumed to be negligible. 
However, the presence of magnitude errors can have a signif- 
icant effect on the shape of the LF. More advanced methods 
for estimating the LF can take this into account (Blanton 
et al. 2003a) but these have not yet been extended to mul- 
tivariate methods. 

Peculiar motions have not been considered in this work. 
However, in order to probe fainter absolute magnitudes, it 
is necessary to reduce the low-redshift limit so that peculiar 
velocities can no longer be neglected. 

A larger sample would allow a more comprehensive cen- 
sus of many galaxy properties. It is important to have good 
galaxy colours, both NIR colours between the four bands 
in the LAS, YJHK, which could reflect various physical 
properties of the galaxies (Eminian et al. 2008), and the 
optical-NIR colours, which would make template fitting pos- 
sible across the nine bands in SDSS and the LAS, which in 
turn would lead to better M/L ratios and stellar masses. 
Neither of these are currently possible, since the Petrosian 
apertures are different in the four LAS bands, and differ- 
ent again to the (single) SDSS Petrosian aperture for each 
galaxy, and since model, e.g. Sersic (1968) profile, magni- 
tudes are not available. Extrapolated total magnitudes such 
as these would also allow a better determination of the bright 
end of the LFs and SMF. 

Sersic (1968) profile fits would also provide structural 
information about the galaxies, as would bulge-to-disc de- 
composition (Driver et al. 2007a,b). 

The limits in radius and surface-brightness also require 
attention. A deficit in low-surface brightness galaxies can be 
accounted for by assuming a certain distribution, e.g. Gaus- 
sian in surface brightness in bins of absolute magnitude, 
and then extrapolating. But it is likely that a lower surface 
brightness limit could be reached using source extraction op- 
timized for low-redshift extended sources, for example, using 
elliptical apertures. Forthcoming deep near-infrared surveys, 
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such as the VISTA (Visible and Infrared Survey Telescope 
for Astronomy) Kilo-degree Infrared Galaxy (VIKING) sur- 
vey, are anticipated to provide measurements to a fainter 
surface brightness limit. 

The colour-dependent incompleteness at low luminosity 
is a result of having a survey limited in both r and K. This 
could be improved by supplementing SDSS redshifts with 
spectroscopy of galaxies within some if-band completeness 
limit. Redshifts from the recently started GAM A 4 spectro- 
scopic survey could be used to probe regions of the param- 
eter space not currently sampled. 



7 CONCLUSIONS 

We have presented the first statistical analysis of galaxies 
from the UKIDSS LAS. The l/Umax space density estimator 
has been used in a four-dimensional form to produce results 
for the K- and r-band LF and the SMF consistent with 
previous findings, with a form similar to a Schechter function 
with an almost fiat faint-end slope. 

We have presented the first iC-band BBD in if-band 
absolute magnitude and effective surface brightness. This 
shows similar trends to the optical BBD: a correlation be- 
tween luminosity and surface brightness, with a broadening 
of the surface brightness distribution at low luminosity and 
a flattening of the luminosity-surface brightness relation at 
high luminosity. 

The multiple limits on the survey have been taken into 
account. For example, limits in if -band and r-band magni- 
tude, if -band Petrosian radius and K -band surface bright- 
ness have been used to estimate the volume within which 
each galaxy would have been visible. 

When the sample is subdivided according to colour, we 
have found a clear distinction between two populations, one 
of red, high-luminosity and high-surface-brightness galax- 
ies, and the other of blue, low-luminosity and low-surface- 
brightness galaxies. This agrees well with previous investi- 
gations of the bimodality of the galaxy population. 

We have considered various sources of incompleteness 
of the data and find that in order to obtain better results 
for low luminosity galaxies with UKIDSS, various selection 
effects need further attention. A stable multivariate estima- 
tor for the space density needs to be developed that takes 
consideration of large-scale structure. Incompleteness at low 
surface brightness also needs to be addressed, perhaps using 
source extraction with elliptical apertures, and the colour- 
dependent selection effects could be addressed by a ff-band 
limited redshift survey to extend the SDSS spectroscopic 
sample to a greater depth in K. Physical properties of galax- 
ies could be investigated further with good NIR and optical- 
NIR colours and/or extrapolated model magnitudes. 



ACKNOWLEDGMENTS 

We thank Seb Oliver for useful discussions, Simon Driver 
and Celine Eminian for comments on an earlier draft of this 
paper, Nigel Hambly and others at the Edinburgh WFAU 

4 Galaxy And Mass Assembly, http://www.eso.org/~jliske/ 
gama/ 



for assistance with the UKIDSS data and Serena Bertone for 
providing model luminosity and stellar mass functions. We 
are grateful to Nigel Metcalfe for providing number count 
data from the literature on his web page. We thank the 
anonymous referee for helpful and constructive comments. 
AJS has been supported by an STFC grant. 

This work is based largely on data obtained as part of 
the UKIRT Infrared Deep Sky Survey. 

Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating In- 
stitutions, the National Science Foundation, the U.S. De- 
partment of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council 
for England. The SDSS Web Site is http : //www . sdss . org/. 

The SDSS is managed by the Astrophysical Research 
Consortium for the Participating Institutions. The Partic- 
ipating Institutions are the American Museum of Natu- 
ral History, Astrophysical Institute Potsdam, University of 
Basel, Cambridge University, Case Western Reserve Uni- 
versity, University of Chicago, Drexel University, Fermilab, 
the Institute for Advanced Study, the Japan Participation 
Group, Johns Hopkins University, the Joint Institute for 
Nuclear Astrophysics, the Kavli Institute for Particle As- 
trophysics and Cosmology, the Korean Scientist Group, the 
Chinese Academy of Sciences (LAMOST), Los Alamos Na- 
tional Laboratory, the Max-Planck-Institute for Astronomy 
(MPIA), the Max-Planck-Institute for Astrophysics (MPA), 
New Mexico State University, Ohio State University, Uni- 
versity of Pittsburgh, University of Portsmouth, Princeton 
University, the United States Naval Observatory, and the 
University of Washington. 

This research has made use of NASA's Astrophysics 
Data System. 



REFERENCES 

Adelman-McCarthy J. K., et al., 2007, ApJS, 172, 634 
Andreon S., 2002, A&A, 382, 495 

Baldry I. K., Glazebrook K., Driver S. P., 2008, MNRAS, 
388, 945 

Ball N. M., Loveday J., Brunner R. J., Baldry I. K., 

Brinkmann J., 2006, MNRAS, 373, 845 
Bell E. F., de Jong R. S., 2001, ApJ, 550, 212 
Bell E. F., Mcintosh D. H., Katz N., Weinberg M. D., 2003, 

ApJS, 149, 289 
Bertone S., De Lucia G., Thomas P. A., 2007, MNRAS, 

379, 1143 

Blanton M. R., et al., 2001, AJ, 121, 2358 
— , 2003a, ApJ, 592, 819 

Blanton M. R., Lin H., Lupton R. H., Maley F. M., Young 
N., Zehavi I., Loveday J., 2003b, AJ, 125, 2276 

Blanton M. R., Lupton R. H., Schlegel D. J., Strauss M. A., 
Brinkmann J., Fukugita M., Loveday J., 2005, ApJ, 631, 
208 

Blanton M. R., Roweis S., 2007, AJ, 133, 734 

Bower R. G., Benson A. J., Malbon R., Helly J. C, Frenk 

C. S, Baugh C. M., Cole S., Lacey C. G., 2006, MNRAS, 

370, 645 

Bruzual G., 2007, in Vazdekis A., Peletier R. F., eds, Proc. 



K-band galaxies from UKIDSS LAS 17 



IAU Symp. 241, Stellar Populations as Building Blocks ol 
Galaxies. Cambridge Univ. Press, Cambridge, p. 125 

Bruzual C, Chariot S., 2003, MNRAS, 344, 1000 

Casali M., et al., 2007, A&A, 467, 777 

Chabrier C, 2003, PASP, 115, 763 

Choloniewski J., 1985, MNRAS, 214, 197 

Cirasuolo M., et al., 2007, MNRAS, 380, 585 

Cole S., et al., 2001, MNRAS, 326, 255 

Cross N., Driver S. P., 2002, MNRAS, 329, 579 

Cross N., et al., 2001, MNRAS, 324, 825 

Croton D. J., et al., 2006, MNRAS, 365, 11 

De Lucia C, Blaizot J., 2007, MNRAS, 375, 2 

Driver S. P., 1999, ApJ, 526, L69 

Driver S. P., et al., 2006, MNRAS, 368, 414 

Driver S. P., Allen P. D., Liske J., Graham A. W., 2007a, 
ApJ, 657, L85 

Driver S. P., Liske J., Cross N. J. G., De Propris R., Allen 

P. D., 2005, MNRAS, 360, 81 
Driver S. P., Popescu C. C, Tuffs R. J., Liske J., Graham 

A. W., Allen P. D., de Propris R., 2007b, MNRAS, 379, 

1022 

Dye S., et al., 2006, MNRAS, 372, 1227 
Efstathiou G., Ellis R. S., Peterson B. A., 1988, MNRAS, 
232, 431 

Eke V. R., Baugh C. M., Cole S., Frenk C. S., King H. M., 

Peacock J. A., 2005, MNRAS, 362, 1233 
Elston R. J., et al., 2006, ApJ, 639, 816 
Eminian C, Kauffmann G., Chariot S., Wild V., Bruzual 

G., Rettura A., Loveday J., 2008, MNRAS, 384, 930 
Gardner J. P., Sharpies R. M., Carrasco B. E., Frenk C. S., 

1996, MNRAS, 282, LI 

Gardner J. P., Sharpies R. M., Frenk C. S., Carrasco B. E., 

1997, ApJ, 480, L99 

Glazebrook K., Peacock J. A., Miller L., Collins C. A., 

1995, MNRAS, 275, 169 
Gott J. R. I., Juric M., Schlegel D., Hoyle F., Vogeley M., 

Tegmark M., Bahcall N., Brinkmann J., 2005, ApJ, 624, 

463 

Graham A. W., Driver S. P., Petrosian V., Conselice C. J., 
Bershady M. A., Crawford S. M., Goto T., 2005, AJ, 130, 
1535 

Hambly N. C, et al., 2008, MNRAS, 384, 637 

Hewett PC., Warren S. J., Leggett S. K., Hodgkin S. T., 

2006, MNRAS, 367, 454 
Huang J.-S., Glazebrook K., Cowie L. L., Tinney C, 2003, 

ApJ, 584, 203 
Huang J.-S., et al., 2001, A&A, 368, 787 
Jarrett T. H., Chester T., Cutri R., Schneider S., Skrutskie 

M., Huchra J. P., 2000, AJ, 119, 2498 
Jones D. H., Peterson B. A., Colless M., Saunders W., 2006, 

MNRAS, 369, 25 
Kochanek C. S., et al., 2001, ApJ, 560, 566 
Lawrence A., et al., 2007, MNRAS, 379, 1599 
Loveday J., 2000, MNRAS, 312, 557 

Lupton R., 1993, Statistics in theory and practice. Prince- 
ton, N.J.: Princeton University Press 

Mailer A. H., Berlind A. A., Blanton M. R., Hogg D. W., 
2009, ApJ, 691, 394 

Mannucci F., Basile F., Poggianti B. M., Cimatti A., Daddi 
E., Pozzetti L., Vanzi L., 2001, MNRAS, 326, 745 

Maraston C, 1998, MNRAS, 300, 872 



Mobasher B., Sharpies R. M., Ellis R. S., 1993, MNRAS, 
263, 560 

Panter B., Heavens A. F., Jimenez R., 2004, MNRAS, 355, 
764 

Petrosian V., 1976, ApJ, 209, LI 

Sandage A., Tammann G. A., Yahil A., 1979, ApJ, 232, 
352 

Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 
525 

Schmidt M., 1968, ApJ, 151, 393 

Sersic J. L., 1968, Atlas de Galaxias Australes. Observato- 

rio Astronomico, Cordoba 
Strauss M. A., et al., 2002, AJ, 124, 1810 
Szokoly G. P., Subbarao M. U., Connolly A. J., Mobasher 

B., 1998, ApJ, 492, 452 



